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ABSTRACT 

Modeling  woven  fabrics  is  challenging  due  to  the  need 
to  capture  the  fabric  response  at  both  the  macroscopic 
scale  of  the  fabric  and  the  “mesoscale”  of  the  yarns  that 
compose  the  weave.  Here,  we  utilize  a  continuum- level 
modeling  technique  for  woven  fabrics  (King  et  ah,  2005) 
which  captures  the  macroscopic  response  of  the  fabric 
as  well  as  the  evolution  of  its  mesostructure  with  defor¬ 
mation  in  a  single  efficient  modeling  step.  A  membrane 
element  formulation  of  the  model  is  implemented  into  a 
nonlinear  finite  element  code  (“Adlib”,  Deiterding  et  ah, 
2006)  with  a  fully  explicit  solver  and  hard  contact.  In 
addition  to  displacement  degrees  of  freedom,  degrees  of 
freedom  representing  the  crimp  amplitude  of  each  yarn 
family  are  defined  in  order  to  solve  for  the  configuration 
of  the  fabric  unit  cell  at  each  step  in  time.  The  accuracy 
of  the  model  is  verified  by  comparing  the  simulation  re¬ 
sults  with  the  experiments  of  Weerasooriya  et  al.  (2008). 

1  INTRODUCTION 

Starting  in  the  1950’s  with  vests  constructed  from 
ballistic  nylon,  fabrics  woven  from  high  strength  fibers 
have  been  crucial  to  the  safety  of  military  personnel. 
The  most  well-known  application  of  woven  fabrics  is 
lightweight,  flexible  body  armor  for  protection  against 
bullets  and  fragmenting  projectiles  from  sources  such  as 
grenades,  mortars,  mines,  and  IED’s.  Woven  fabrics  are 
also  used  to  reinforce  polymers  or  ceramics  in  helmets 
and  armor  panels  for  protection  against  high-velocity 
projectiles.  Additional  common  uses  include  restraint 
systems  such  as  seat  belts  and  head  restraints,  deploy¬ 
able  fabric  structures  including  air  bags  and  air  beams, 
as  well  as  other  forms  of  ballistic  protection  such  as  pilot 


seats.  Most  recently,  fabrics  are  being  utilized  to  increase 
the  portability  and  durability  of  personal  electronic  de¬ 
vices. 

Despite  decades  of  research,  the  design  of  woven  fabric 
systems  is  still  largely  based  on  empirical  methods  and 
experience.  Modeling  woven  fabrics  is  challenging  due 
to  the  need  to  capture  the  fabric  response  at  both  the 
macroscopic  scale  of  the  fabric  and  the  “mesoscale”  of 
the  yarns  that  compose  the  weave.  In  fact,  there  is  no 
widely  accepted  computational  model  for  fabric  defor¬ 
mation  which  is  both  computationally  efficient  and  able 
to  represent  the  behavior  and  evolving  structure  of  the 
fabric  at  the  level  of  the  individual  yarns.  In  the  litera¬ 
ture,  either  the  fabric  has  been  modeled  at  the  level  of 
the  mesoscale,  with  detailed  discretization  of  each  fiber 
in  the  weave,  or  the  fabric  has  been  modeled  at  the  con¬ 
tinuum  level,  with  phenomenological  approaches,  typ¬ 
ically  treating  the  fabric  as  an  anisotropic  continuum 
with  two  preferred  material  directions.  In  the  first  ap¬ 
proach,  the  discretization  length  scale  is  determined  by 
the  dimensions  of  the  fibers  and  weave,  and  thus  real¬ 
sized  or  multi-ply  simulations  are  computationally  pro¬ 
hibitive.  The  second  approach  cannot  be  used  to  design 
fabrics  based  on  the  properties  of  the  fibers  and  weave, 
and  it  is  unsuitable  for  applications  where  the  macro¬ 
scopic  response  depends  on  the  fabric’s  unique  under¬ 
lying  mesostructure  and  its  specific  deformation  mech¬ 
anisms,  such  as  fiber  failure,  crimp  interchange,  shear 
locking,  or  relative  slip  between  the  two  fiber  families.  As 
such,  without  suitable  models,  the  analysis  and  improve¬ 
ment  of  existing  products  and  the  development  of  new 
applications  for  woven  fabrics  is  currently  overly  costly 
and  time  consuming. 

Here,  we  utilize  a  continuum-level  modeling  technique 
for  woven  fabrics  (King  et  ah,  2005)  which  captures  the 
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macroscopic  response  of  the  fabric  as  well  as  the  evo¬ 
lution  of  its  mesostructure  with  deformation  in  a  single 
efficient  modeling  step.  An  appropriate  unit  cell  approx¬ 
imation,  based  on  a  pin-joined  truss  geometry,  represents 
the  weave  architecture.  The  configuration  of  the  unit  cell 
is  determined  by  the  macroscopic  deformation  gradient 
and  the  satisfaction  of  local  equilibrium.  The  forces  on 
the  unit  cell,  and  in  turn  the  corresponding  continuum- 
level  stresses,  are  calculated  from  physically  motivated 
constitutive  relations. 

The  fabric  model  is  used  to  simulate  projectile  impact 
on  single  plies  of  plain- weave  Kevlar.  Direct  comparisons 
are  made  with  experimental  data  obtained  from  digital 
image  correlation  (DIC). 

2  FABRIC  MODEL 

The  woven  fabric  model  of  King  et  al.  (2005)  has  been 
modified  and  extended  in  order  to  function  efficiently 
within  an  explicitly  integrated  finite  element  framework 
(“Adlib”,  Deiterding  et  al.,  2006). 

2.1  Constitutive  Model 


Figure  1:  Geometry  of  fabric  unit  cell  (King  et  al.,  2006). 

The  fundamental  basis  of  the  model  is  the  assumed 
unit  cell  geometry  of  the  fabric  shown  in  Figure  1.  The 
unit  cell  approximates  the  mesoscale  mechanisms  of  de¬ 
formation  that  underlie  the  macroscopic  response  of  the 
fabric.  The  continuum  response  at  each  material  point 
along  the  fabric  surface  is  obtained  through  homogeniza¬ 
tion  of  the  unit  cell  response,  which  is  in  turn  calculated 
through  the  analytical  model  described  below.  The  warp 
(k  =  1)  and  weft  (k  =  2)  yarns  are  approximated  as  pin- 
joined,  elastic  trusses  with  the  following  properties: 
pk  yarn  quarter- wavelength 
Lk  yarn  half-length  or  simply  “length” 


Figure  2:  Geometry  of  fabric  unit  cell  with  locking  spars 
(King  et  al.,  2006). 


f3k  yarn  crimp  angle 
Ak  yarn  crimp  amplitude 
0  in-plane  angle  between  yarns 

The  torsional  bending  spring  provides  the  resistance  to 
bending  of  the  yarns  while  the  cross-over  spring  provides, 
at  the  cross-over  points,  torsional  resistance  to  in-plane 
shear  as  well  as  compressive  resistance  to  compaction 
of  the  yarns.  Figure  2  illustrates  the  locking  trusses  of 
length  dk  and  angle  ak  from  the  plane  of  the  fabric  which 
resist  both  shear  and  cross  locking  of  the  yarns.  Consti¬ 
tutive  relations  describe  each  mode  of  deformation  of  the 
yarns.  They  are  (with  superscript  “0”  indicating  an  ini¬ 
tial,  undeformed  value): 


1. 

2. 

3. 


4. 


5. 


extension 

yarn  tension,  T k  =  kk(Lk  —  Lk) 
bending 

bending  moment,  Mbk  =  kbk(f3k  Pk) 


compaction 

force  in  cross-over  spring,  F\  =  Ah(eaI  —  1),  where 

I  =  —  AAi  —  AA2  is  the  interference  of  the  yarns  at 

the  cross-over  points. 

yarn  locking 

force  in  locking  truss, 


0  dk  —  dk  <  0 

Kd(d°k-dky  d°k-dk>  o 


in-plane  shear 

moment,  Ms  =  As7e,  where  the  total  relative  yarn 
rotation,  7  =  6°  —  6,  is  separated  into  elastic  and 
dissipative  components,  7  =  7e  +  qf,  which  evolve 


according  to  7f  =  70 


fcfc,  fcbfej  Ah,  a,  Ad,  c,  7o,  and  Mo  are  constants  deter¬ 
mined  experimentally  from  simple  experiments  on  single 
yarns  and  strips  of  fabric. 

The  deformed  geometry  of  the  unit  cell  in  three  dimen¬ 
sions  is  fully  determined  by  the  macroscopic  deformation 
gradient,  F  =  Grad  x(X,  £),  and  through-thickness  equi¬ 
librium.  F  specifies  the  orientation  and  stretch  of  the 
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yarn  vectors,  gi  and  g2,  everywhere  in  the  fabric: 


2.2  Finite  Element  Model 


gi  =  Fg?  and  g2  =  Fg^.  (1) 


Pi,  P2 1  and  0  are  calculated  from  gi  and  g2-  The  out- 
of-plane  state  of  the  fabric  (crimp  amplitude)  is  derived 
from  a  force  balance  at  the  cross-over  spring.  At  the 
end  of  each  cross-over  spring,  for  static  equilibrium,  the 
force  in  the  cross-over  spring,  Fj,  must  balance  the  net 
through-thickness  component  of  the  force  in  the  yarn 
trusses,  F^: 

Fi  =  Fk,  k  =  l,2.  (2) 

Fj  is  given  above,  and  Fj^1  depends  on  the  moments  in  the 
bending  springs  and  the  forces  in  the  yarns  and  locking 
spars: 


Fk 


27i,.  sin/3/,.  —  2FL;siiia:/ 

kbkjPk  ~  P®) 


1  + 


Pi\cos9\cos/3k\ 


Lk 


+2 


Lk 


cospk, 


(3) 


where  /  =  1,  2  (/  ^  k). 

Finally,  as  illustrated  schematically  in  Figure  3,  in  or¬ 
der  to  homogenize  the  mesoscopic,  yarn-level  forces  to 
continuum- level  stress  measures,  the  stress  on  the  unit 
cell  must  be  computed.  The  first  Piola-Kirchoff  stress, 
P  =  P(F,  AAi,  AA2),  is  therefore  calculated  by  equili¬ 
brating  the  external  stress  which  must  act  on  the  loaded 
unit  cell  with  the  in-plane  forces  that  the  yarns  exert  on 
the  unit  cell. 


Continuum  Approximation 


The  fabric  model  has  been  integrated  with  the  total  La- 
grangian  framework  of  an  explicit  finite  element  code 
(“Adlib”,  Deiterding  et  ah,  2006),  in  order  to  obtain  a 
versatile  tool  for  simulations  of  ballistic  impacts  on  wo¬ 
ven  fabrics  (WovenShield,  Parasim  Inc.).  The  degrees  of 
freedom,  d,  are  the  three  components  of  macroscopic  dis¬ 
placement,  Ui:  and  the  changes  in  yarn  crimp  amplitude, 
AAi  and  AA2.  The  degrees  of  freedom  are  explicitly  in¬ 
tegrated  forward  in  time  via  the  Newmark  method.  At 
£n+1, 


and 


At2  ■■ 

dn+i  =  dn  -f-  Atdn  —  d„, 

(4) 

dn+i  =  M-1  (f T+l+fntl), 

(5) 

dn+i  =  dn  -| - ^dn  +  dn+i^  . 

(6) 

The  mass  matrix,  M,  is  diagonal  with  the  mass  appor¬ 
tioned  equally  throughout  the  mesh. 

The  majority  of  resistance  to  ballistic  impact  loads  is 
generated  from  in-plane  deformation  of  the  fabric  and 
out-of-plane  inertial  effects.  Out-of-plane  bending  and 
shear  provide  negligible  contributions  to  the  ballistic  per¬ 
formance.  As  such,  it  was  deemed  sufficient  to  model 
the  fabric  as  a  membrane.  6-node  triangular  and  9-node 
rectangular  elements  were  used. 

The  in-plane  position  of  a  material  point  in  the  refer¬ 
ence  configuration  is  given,  as  a  function  of  the  natural 
coordinates,  (s  =  1,  3  for  triangular  elements;  s  =  1,2 
for  rectangular  elements)  by  the  shape  functions  at  each 
node,  a: 


Xj(Zs)  =  XjaNa(Zs\  J  =  1,2.  (7) 

Correspondingly,  the  three-dimensional  position  of  a  ma¬ 
terial  point  in  the  current  configuration  is  given  by 

=  XiaNa(£s),  i=  1,3.  (8) 

The  relevant  components  of  the  deformation  gradient  at 
a  point  within  the  element  are  therefore  defined  by 

FiJ  =  ~FV~  =  N*,JXia-  (9) 

Provided  that  the  fabric  initially  lies  in  the  1-2  plane,  F 
then  determines  the  orientation  and  stretch  of  the  yarn 
vectors,  gi  and  g2,  at  each  integration  point  within  the 
element.  The  internal  nodal  forces  corresponding  to  the 
displacement  degrees  of  freedom  are  then  calculated  by 
integrating  numerically  over  the  element: 


Figure  3:  Modeling  fabric  as  a  homogeneous  anisotropic 
continuum  with  current  yarn  directions  gi  and  g2  (King 
et  ah,  2006). 


(/Lnt)e  =  —  fn0e  Fij(F)NatjdVo.  (10) 

By  virtual  work  considerations,  multiplying  Equa¬ 
tion  2  by  virtual  changes  in  amplitude  Ai^  and  inte¬ 
grating  over  the  element  provides  the  statement  of  equi¬ 
librium  (excluding  inertia  and  external  forces)  in  weak 
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form: 


AA, 


FjjA^AA^)  -  FfcN( F,  AAUAA2) 

4 P1P2 


dS0  =  0. 


(11) 

With  A Ak(^s)  =  AAkaNa(£,s),  the  internal  nodal  forces 
corresponding  to  the  amplitude  degrees  of  freedom  be¬ 
come 


(  ^int\e  _ 

\J ka  )  — 


F1(AA1,AA2)  -  FfeN( F,  AAi,AA2) 


W  W2 

The  total  mass  of  the  element  is  given  by 

3 

.  Wq' 

2 


1  3 

M  =  pftf-W"YJJq> 


NadS0. 

(12) 

(13) 


q=  1 


where  pf  and  R  are  the  density  and  thickness,  respec¬ 
tively,  of  the  fabric  and  Jq  and  Wq  are  the  Jacobian  de¬ 
terminant  and  weight,  respectively,  at  integration  point 
q.  In  order  to  calculate  the  mass  associated  with  the 
amplitude  degrees  of  freedom,  we  take  L\  «  =  T, 

(3\  ~  =  (3  s  0  and  m\  ~  m2  =  m  =  M/4.  Moment 

equilibrium  of  a  single  half-length  of  yarn  about  its  point 
of  intersection  with  the  unit  cell  boundary  gives 


the  experimental  results  of  King  et  al.  (2005)  for  uniaxial 
tension  of  2.54  cm  x  25.4  cm  strips  of  fabric.  The  fab¬ 
ric  tested  by  King  et  ah,  commonly  known  as  S706,  was 
600  denier,  plain-weave  Kevlar  KM2  with  34  threads  per 
inch  in  both  the  warp  and  weft  directions.  All  parame¬ 
ters  used  in  the  model  were  either  measured  by  King  et 
al.  or  fit  to  the  experiments  of  King  et  al.  640  6-node 
planar  triangular  elements  were  used  to  model  the  strip. 
Excellent  agreement  is  seen  for  loading  in  both  the  warp 
and  weft  directions. 


FI  «  1(3 


—p. 

12  M 


(14) 


With  the  through-thickness  acceleration  of  the  free  end 
approximated  as  Ak  ~  (31,  the  through-thickness  equilib¬ 
rium  equation  becomes 


F  = 


M  •• 

— Ak  • 
12 


(15) 


Thus,  when  the  mass  of  the  other  yarn  half-length  is  in¬ 
cluded,  we  find  that  the  lumped  masses  associated  with 
the  amplitude  degrees  of  freedom  are  related  to  the  cor¬ 
responding  masses  associate  with  displacement  degrees 


of  freedom  by  a  factor  of  -. 
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Figure  5:  Simulations  vs.  experiments  (King  et  al.,  2005) 
for  strips  of  Kevlar  S706  under  uniaxial  tension  in  warp 
and  weft  directions. 


Figure  6  plots  simulation  and  experimental  results  for 
bias-extension.  In  this  case,  a  3.5  cm  x  9.5  cm  strip  of  fab¬ 
ric,  with  the  yarns  aligned  at  45°  to  the  horizontal,  was 
loaded  in  uniaxial  tension.  300  6-node  planar  triangular 
elements  were  used  to  model  the  strip.  The  remarkable 
agreement  between  the  simulation  and  experiments  pro¬ 
vides  a  validation  of  the  ability  of  the  model  to  simulate 
complex  deformations. 

3  EXPERIMENTAL  PROCEDURE 


Figure  4:  Calculation  of  lumped  mass  corresponding  to 
amplitude  degrees  of  freedom. 


2.3  Quasi-Static  Model  Verification 

In  order  to  verify  the  accuracy  of  this  new  formulation  of 
the  fabric  model,  several  standard  test  cases  were  simu¬ 
lated.  In  Figure  5,  simulation  results  are  compared  with 


Weerasooriya  et  al.  (2008)  conducted  projectile  im¬ 
pact  experiments  on  single  plies  of  the  Kevlar  S706. 
A  gas  gun  was  used  to  launch  d  =  12.7  mm 

hemispherically-nosed  steel  projectiles  of  mass  104  g  at 
a  striking  velocity  of  22  m/s.  As  shown  in  Figures  7 
and  8,  the  approximately  40  cm  square  fabric  specimens 
were  stapled  to  30.5  cm  square  “picture- frames”  with  a 
pneumatic  staple  gun.  Frames  with  a  25.4  mm  square 
cut-out  and  a  R  =  12.7  cm  circular  cut-out  were  used. 
The  frames  were  clamped  between  two  aluminum  plates 
of  similar  dimensions.  The  surface  of  the  aluminum  plate 
in  contact  with  the  fabric  was  roughened  in  order  to  limit 
slipping  of  the  fabric.  The  specimens  were  first  painted 
with  white  spray  paint,  and  then  a  speckle  pattern  was 
applied  with  black  spray  paint. 
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Figure  6:  Simulation  vs.  experiments  (King  et  al.,  2005) 
for  Kevlar  S706  bias-extension  test. 


Figure  9  illustrates  the  experimental  set-up.  The 
launch  velocity  of  the  projectile  was  measured  at  the 
outlet  of  the  gun  barrel  with  laser  beam  detectors. 
The  detectors  also  triggered  the  two  Photron  APX-RS 
high  speed  digital  cameras  to  start  recording  images. 
256  x  256  pixel  images  of  the  back  surface  of  the  speci¬ 
men  were  acquired  at  a  rate  of  30  kHz.  Commercial  digi¬ 
tal  image  correlation  software  from  Correlated  Solutions, 
Inc.  was  then  used  to  measure  the  full-field  in-plane  (u, 
v )  and  transverse  (w)  displacements  of  the  fabric. 

4  RESULTS 

The  fabric  model  was  used  to  simulate  the  experi¬ 
ments  of  Weerasooriya  et  al.  2912  6-node  triangular 
membrane  elements  and  900  9-node  rectangular  mem¬ 
brane  elements,  respectively,  were  used  in  the  circular 


Figure  7:  Front  view  of  Kevlar  fabric  specimen  stapled 
to  frame  with  square  hole. 


Figure  8:  Back  view  of  Kevlar  fabric  specimen  stapled 
to  frame  with  round  hole. 


and  square  cases.  Contact  between  the  projectile  and 
fabric  was  simulated  as  hard  and  without  friction.  While 
great  care  was  taken  in  the  experimental  tests  to  mini¬ 
mize  slip  of  the  fabric  from  the  clamps,  as  with  nearly 
all  experiments  of  this  type,  some  slip  did  occur.  For 
example,  the  DIC  results  for  in-plane  u-displacement  in 
Figure  10  show  that,  during  the  R  =  12.7  cm  experiment, 
at  least  4  mm  of  slip  occurred  from  the  clamps.  On  av¬ 
erage,  3-5  mm  of  slip  is  apparent  for  both  the  square  and 
circular  tests.  Therefore,  in  the  simulations  the  edges  of 
the  specimen  were  not  considered  perfectly  fixed  at  the 
clamps.  Rather,  a  distributed  frictional  force  consistent 
with  the  observed  slip  of  the  fabric  was  imposed  on  the 
perimeter  of  the  fabric. 

In  Figures  11-14,  the  transverse  displacements  along 
lines  in  the  warp  and  weft  directions  through  the  point  of 
impact  of  the  two  fabrics  are  plotted  for  the  experiments 
at  four  instances  in  time,  with  the  final  time  correspond- 


SDecimen  back  surface 


Figure  9:  Schematic  of  experimental  set-up  for  measur¬ 
ing  full- field  displacements  during  a  ballistic  impact  test. 
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Figure  10:  Contours  of  in-plane  u- displacement  during 
R  =  12.7  cm  experiment  at  t  —  2.0  ms. 


ing  to  the  maximum  displacement  of  the  fabric.  During 
the  experiment,  the  fabric  deforms  almost  identically  in 
the  warp  and  weft  directions.  The  isotropy  of  the  ex¬ 
perimental  results  is  an  indication  that  the  specimen  is 
subject  to  a  small  prestretch.  In  fact,  it  is  nearly  impos¬ 
sible  to  clamp  flat  a  fabric  specimen  without  applying 
at  least  a  small  degree  of  preload.  In  the  simulations, 
the  effects  of  preload  are  accounted  for  by  specifying  an 
initial  deformation  gradient,  Fpre,  and  initial  changes 
in  yarn  amplitude,  AA^  and  A  A®.  These  parameters 
were  determined  by  subjecting  a  single  4-node,  reduced 
integration  finite  element  to  increasing  levels  of  biaxial 
tension.  Single  element  simulations  were  conducted  with 
both  equal  force  and  equal  strain  in  the  warp  and  weft 
directions.  A  parametric  study  with  varying  levels  of 
Fpre  and  A  A®  was  conducted  for  the  R  =  5  inch  simu¬ 
lations.  It  was  found  that  a  preload  with  equal  force  in 


Figure  12:  Experiment  versus  simulation  of  R  =  12.7  cm 
fabric,  clamped  at  circumference,  for  V{  =  22  m/s:  weft 
yarn  transverse  displacement  through  point  of  impact. 


the  warp  and  weft  directions  provided  a  good  fit  to  the 
experimental  data.  The  preload  was  found  to  be  1.7  N 
per  yarn,  and  the  corresponding  initial  changes  in  ampli¬ 
tude  were  AA^  =  —0.013  mm  and  AAl]  =  —0.031  mm. 
The  results  of  the  simulations  with  both  fabric  slip  at 
the  clamps  and  preload  compared  with  the  experiments 
are  also  plotted  in  Figures  11-14.  Exceptional  agreement 
between  experiment  and  simulation  is  exhibited  for  both 
geometries  of  fabric  specimen. 

Figure  15  plots  the  history  of  displacement  of  the  fab¬ 
ric  at  the  point  of  impact  for  both  the  experiments  and 
simulations.  Overall,  the  simulations  compare  very  well 


Figure  13:  Experiment  versus  simulation  of  25.4  cm 
square  fabric,  clamped  on  four  sides,  for  V[  =  22  m/s: 
warp  yarn  transverse  displacement  through  point  of  im¬ 
pact. 


Figure  11:  Experiment  versus  simulation  of  R  =  12.7  cm 
fabric,  clamped  at  circumference,  for  v\  =  22  m/s:  warp 
yarn  transverse  displacement  through  point  of  impact. 
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5  CONCLUSIONS 


Figure  14:  Experiment  versus  simulation  of  25.4  cm 
square  fabric,  clamped  on  four  sides,  for  v\  =  22  m/s: 
weft  yarn  transverse  displacement  through  point  of  im¬ 
pact. 


with  the  experiments,  even  during  the  rebound  phase  of 
the  deformation.  Figure  16  shows  the  deformed  shape 
and  contours  of  transverse  displacement  for  the  experi¬ 
ments  and  simulations  with  R  =  12.7  cm.  As  expected 
from  the  previous  results,  there  is  an  impressive  cor¬ 
relation  between  the  two.  Furthermore,  we  note  the 
“tent-like”  deformation  patterns  of  both  the  experiment 
and  simulation.  Observed  by  numerous  experimentalists, 
this  mode  of  deformation  is  due  to  the  highly  anisotropic 
properties  of  the  fabric. 


Figure  15:  Experiment  and  simulation  transverse  dis¬ 
placement  at  point  of  impact  for  =  22  m/s:  R  = 
12.7  cm  and  25.4  mm  square. 


The  extension  of  the  fabric  model  to  an  explicit  finite 
element  formulation  with  the  deformed  configuration  of 
the  unit  cell  determined  solely  from  degrees  of  freedom 
driven  by  inertia  has  been  a  success.  The  model  has 
proven  to  simulate  projectile  impact  events  accurately 
without  excessive  computational  cost. 

The  model  has  many  potential  uses.  Its  computational 
efficiency  enables  the  simulation  of  large  and  complex 
(including  multi-ply)  structures.  Although  a  continuum- 
level  model,  it  also  simultaneously  monitors  the  deforma¬ 
tion  of  the  fabric  on  the  scale  of  individual  yarns.  Knowl¬ 
edge  of  local  quantities  such  as  yarn  tension  and  yarn 
wavelength  provides  the  ability  to  predict  phenomenon 
such  as  penetration  due  to  yarn  breakage  (e.g.  bullet, 
knife  threat)  or  yarn  separation  (puncture  threat). 

The  framework  of  the  model  is  versatile  and  may  be 
applied  to  many  existing  or  proposed  materials.  Other 
plain-weave  Kevlar  fabrics  may  be  modeled  by  simply 
adjusting  the  purely  physical  material  parameters.  Fab¬ 
rics  of  a  different  weave  pattern  require  only  a  new  unit 
cell  geometry.  Viscoelastic  yarn  deformation  may  be  in¬ 
cluded  in  the  model  by  modifying  the  tensile  constitutive 
response  of  the  yarns.  Furthermore,  its  continuum  na¬ 
ture  will  allow  the  model  to  be  integrated  with  other 
material  models  for  the  simulation  of  fabric  composites 
such  as  those  used  in  helmets  and  armor  panels. 

ACKNOWLEDGEMENTS 

This  research  was  supported  by  the  U.S.  Army  Re¬ 
search  Office  through  the  MIT  Institute  for  Soldier 
Nanotechnologies  under  contract  number  W911NF-07- 
D-0004.  The  authors  thank  Professor  R.  Radovitzky  for 
providing  the  explicit  finite  element  framework,  Adlib, 
in  which  we  implemented  the  fabric  model. 

REFERENCES 

Deiterding  R.,  Radovitzky  R.,  Mauch  S.P.,  Noels  L., 
Cummings  J.C.,  and  Meiron  D.I.,  2006:  A  Virtual 
Test  Facility  for  the  Efficient  Simulation  of  Solid  Ma¬ 
terial  Response  Under  Strong  Shock  and  Detonation 
Wave  Loading,  Eng.  Comp.,  22,  325-347. 

King  M.,  2006:  A  Continuum  Constitutive  Model  for 
the  Mechanical  Behavior  of  Woven  Fabrics  Including 
Slip  and  Failure,  PhD  Thesis,  MIT,  Cambridge,  MA. 
King  M.,  Jearanaisilawong  P.,  and  Socrate  S.,  2005:  A 
Continuum  Constitutive  Model  for  the  Mechanical 
Behavior  of  Woven  Fabrics,  Int.  J.  Sol.  Struct .,  42, 
3867-3896. 

Weerasooriya  T.,  Gunnarsson  C.A.,  and  Moy  P.:  Mea¬ 
surement  of  Full-Field  Transient  Deformation  of  the 
Back  Surface  of  a  Kevlar  KM2  Fabric  During  Im¬ 
pact  for  Material  Model  Validation,  Proc.  2008  Int. 
Cong.  Expos.  Exp.  Mech.  Appl.  Mech .,  June  2-5, 
Orlando,  FL. 


7 


Figure  16:  Contours  of  transverse  displacement  for  experiment  (left,  in  mm)  and  simulation  (right,  in  m)  for 
R  =  12.7  cm  fabric  with  v\  =  22.2  m/s  as  a  function  of  time:  (a),  (a')  t  =  0.4ms;  (b),  (b')  t  =  0.9ms;  (c),  (c') 
t  =  1.3ms;  (d),  (d')  t  =  2.0ms. 
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